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We continue the development of a manifestly 4-dimensional, completely covariant, approach to 
transformation optics in linear dielectric materials begun in a previous paper. This approach, which 
generalizes the Plebanski based approach, is systematically applicable for all transformations and all 
general linear materials. Importantly, it enables useful applications such as arbitrary relative motion, 
transformations from arbitrary non-vacuum initial dielectric media, and arbitrary space-times. This 
approach is demonstrated for a resulting material that moves with uniform linear velocity, and in 
particular for a moving cloak. The inverse problem of this covariant approach is shown to generalize 
Gordon's "optical metric". 
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I. INTRODUCTION 

The emerging field of transformation optics, where useful arrangements of man-made "metamaterials" [H-Q are 
designed via transformations of electromagnetic fields, has theoretical roots stretching back nearly a century, to the 
early days of general relativity. The idea that the behavior of light in a gravitational field can be replicated by a 
suitable distribution of refractive media appears to have been first postulated by Eddington . Subsequently, Gordon 
studied the inverse problem, that of representing a refractive medium as a vacuum space-time described by an 
"optical metric". Later, Plebanski found effective constitutive relations for electromagnetic waves propagating 
in vacuum space-times, but, while recognizing the formal equivalence of these equations to those in a macroscopic 
medium in flat space-time, does not exploit this equivalence to actually describe such a medium. This was first done 
by De Felice Q , who used the Plebanski equations to describe the equivalent medium of both a spherically symmetric 
gravitational system and Friedmann-Robertson- Walker space-time. 

More recently, Pendry [8] , pointed out the specific relationship between spatial transformations and material prop- 
erties, and demonstrated how it could be used to create novel devices. Closely related work by Greenleaf et al. Q 
developed a similar concept for electric current flow and applied it to impedance tomography. The initial approaches 
to transformation optics relied on purely spatial transformations [1, [l^, [ll[ . Leonhardt and Philbin generalized this 
to transformations involving both space and time by using the explicit equivalence given by De Felice; thus linking 
transformation optics to differential geometry p^ . Another approach based on field-transforming metamateri- 
als, considers more general transformations in the Fourier domain, but at the expense of the intuitive and appealing 
geometric interpretation. For more recent reviews see Refs. [l^ [la|. 

In examining the details of the Plebanski-De Felice approach, we find some limitations that are addressed here. 
As pointed out by Plebanski himself [6|, the equation now bearing his name is not strictly covariant, because its 
derivation requires a matrix inversion that is not a true tensor operation. One consequence of this, in the context 
of transformation optics, is that magneto-electric coupling terms can not always be simply interpreted as a material 
velocity, as is frequently done. This is because the Plebanski equations can only be identified with stationary media 
or with slowly moving (i.e. nonrelativistic) , isotropic media pi3i |. This prompts the question of whether the approach 



of Ref. |12l | may be generalized to clearly distinguish magneto-electric couplings from material velocity, and allow for 
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specially relativistic corrections. Just such an approach was outlined in Ref. fl^, and was demonstrated to to recover 
several results obtained through other means. 

Here we provide a complete derivation of the approach outlined in Ref. [l6j , further generalizing the result obtained 
there. A physically realistic scenario for transformation optics designed devices is that the device move with arbitrary 
velocity. We demonstrate that the approach described here may be applied to find the material properties of a 
transformation when the resultant material is constrained to move with arbitrary uniform velocity with respect to 
the frame in which the transformations are given. This represents a departure from most previous examples in 
transformation optics, where either the resulting material is stationary or where the velocity is dictated by the 
transformation itself. While a few examples exist of nonrelativistic moving dielectrics in transformation optics for 
special cases [H, [13] , we provide a systematic approach that is widely applicable for any velocity. 

Another limitation of the Plebanski-De Felice approach is that the resulting material must reside in vacuum, 
Minkowski space-time. The approach described here relaxes these conditions, allowing for physically realistic scenarios 
such as transformations in arbitrary non- vacuum initial dielectric media (isj , or in arbitrary space-times - thus 
providing general relativistic corrections for transformations in arbitrary space-times, such as the weakly curved 
space-time near Earth ^1^. Lastly, we show that this covariant approach is consistent with, and generalizes, Gordon's 
optical metric as essentially the inverse problem of transformation optics. 

The paper is organized as follows: In Sec.|lT]we review the completely covariant theory of vacuum electrodynamics 
using modern language; for which the necessary ideas and notation from differential geometry may be found in 
Appendix [B] This review is presented in some detail, because in Sec. IIIII the covariant theory of Sec. |ll]is extended 
to describe electrodynamics in macroscopic linear dielectric materials. This section presents a slight departure from 
the usual description of electrodynamics in dielectric materials in order to clearly distinguish material effects from 
space-time effects. Section HVl describes the concept of transformation optics and presents an interpretation consistent 
with the geometric picture of the preceding sections. The main result for applications in transformation optics is Eq. 
(j30p . Section llV Al examines a particular transformation both when the resulting material is at rest and when it is 
in motion relative to the frame in which the fields have been measured. As expected, it is found that the results for 
the material in motion smoothly recover, in the limit zT — ?> 0, the results for the material at rest. In Sec. |V] we study 
the inverse problem of transformation optics, that of finding an equivalent vacuum space-time starting from an initial 
dielectric, thus generalizing Gordon's optical metric idea. We conclude with Sec. I VII 



II. CLASSICAL ELECTRODYNAMICS IN VACUUM 



The basic elements of covariant electrodynamics needed for transformation optics were presented in Ref. 16], here 
we present a relatively self-contained and more detailed description of covariant electrodynamics in both vacuum 
and linear materials. This development relies on the geometric language and tools of differential geometry, such as 
exterior derivative, wedge product, and the pullback of a tensor. These aspects of differential geometry are described 
in a myriad of excellent sources, such as |20l423l |. The most important of these for our purpose, the pullback map, 
is described in Appendix |BJ The development and notation, in particular the sign convention, follows that of Ref. 
[23 |, while more information, particularly for electrodynamics in materials, can be found in Refs. [1^, [2^. We use the 
Einstein summation convention, indices are lowered (raised) by the metric tensor gap (its inverse 5"^), and the speed 
of light and Newton's constant are set to c = G = 1. 



A. Field Strength Tensor 



In free space, classical electrodynamics is modeled as a principal U{1) fiber bundle over a space-time manifold M 
(which we assume to be equipped with a metric) with connection 1-form A = (frequently called a "gauge field" , 
A is the covariant version of the 4- vector potential). The field strength F = F^i, is the curvature 2-form of the U{1) 
fiber bundle, equal to the exterior derivative of A, 

F = dA ^ F,,, = A,,,^ - A,,, (1) 

where the comma indicates a derivative. The components of F can be represented as a matrix, that in a local 
orthonormal frame (or Minkowski space-time with Cartesian coordinates) have values 
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With this choice of values for the components of F, the 4- force vector on a particle moving with 4- velocity u" and 
charge q is 

r = qg'^^F^uW, (3) 

which is just the Lorentz force. For example, a particle at rest with respect to this system has 4- velocity = (1,0,0,0), 
for which = g(0, E^, Ey, E^), recovering the usual notion that a charged particle at rest feels only the electric part 

of the field. For a particle moving with u" = 7(1, u^, v^,v^), the spatial components of are / — (^E + v x B^, 

while the time component is the change in energy per unit time, or the power. 

The E and B fields are now tightly intertwined, simply representing different components of a single object, F. 
Because the second exterior derivative of any form vanishes and the fact that F = dA, it immediately follows that 

dF = 0. (4) 

This is nothing more than the covariant form of the homogeneous Maxwell equations, and shows that the homogeneous 
equations are simply geometric conditions imposed on the fields. 



B. Field Strength Dual 

Naturally associated to each point of the space-time manifold are four m-dimcnsional vector spaces, where m = 
dim(M). The metric generates a bijection g between the space of fc-forms a'^T*{M) and the space of A:- vectors 
/\^Tp{M). The volume form provides a bijection w between the space of fc-forms and the space of (m — fc)-vectors 
A^'"~'^^Tp(M). The composition of maps, called the Hodge dual, makes the diagram of figure 1 commutative. The 
Hodge dual provides a natural two-form dual to the field strength F^^. In particular, we define the map 



★ : A%*(M) A^T*{M) 



(5) 



as the composition -k = uj o g applied to 2-forms (where A T*(M) is the space of 2-forms) , or in component form 



1 



(6) 



The components of ^F can also be represented as a matrix, that in a local orthonormal frame (or Minkowski space-time 
with Cartesian coordinates) have values 
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The dual nature is now explicit; where we had the decomposition F = Eadx"" A df + Bo6(da;'* A dx**) we now have the 
dual decomposition *F = —BaAx"^ /\&t + Eab{<^x"' A da;**). 



/\{m-k)nip*(^M) -(-^ A('"-*=)Tp(M) 




A%*(M) ^ 



FIG. 1: The metric generates a bijection g between (p) and (°) alternating tensors, while the volume form w provides a bijection 
between (p) and (^^^'1^,) alternating tensors (for any k < m). The composition -k = to o g provides a bijection between ('"^'') 
and (q) alternating tensors and between {^^_^) and (J^) alternating tensors. 
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C. Vacuum Action 



To describe electrodynamics by means of a variational principle requires an action, and to construct an action 
S = J d^x \/\g\C requires a suitable Lagrangian density jC. Differential geometry tells us that we can only integrate 
a fc-form over a /c-dimensional (sub)manifold. There are 3 ways to construct a 4-form from F and ★F that can be 
integrated over the 4-dimensional space-time manifold: F A F, ★F A ★F, and F A ★F (we neglect a possibility such 
as A A -kA, which results in a massive photon described by the Proca equations .27]). Consider the first possibility. 
Since F = dA, 

/ (FAF)= / d(AAF)= / (AAF). (8) 

Jm Jm JdM 

Essentially, this term is a total 4-divergence and contributes nothing to the integral. Writing out the wedge product 
shows that ★F A ^F = — F A F, so the second possibility also contributes nothing to the integral. The last possibility, 
however, gives a contribution 

/ (FA*F)=/ d\^\g{¥,Y)^[ d'x^\{F'^'^F^,). (9) 

JM JM JM 

The component form of this expression - the right hand side of Eq. ([9]) - is the version most commonly encountered in 
the literature. That this appears in component form as F'^'^ Ff^^, should be regarded as a happy coincidence that is an 
artifact of the classical vacuum; do not loose sight of the fact that the fundamental tensors making up the Lagrangian 
density are F and ^F, as this will be generalized shortly. 



D. Excitation Tensor 



The field strength tensor F encodes some information about the fields, namely the electric field strength and the 
magnetic flux. Consider that on the other side of this coin, the magnetic field strength and the electric flux are 
encoded in another tensor G, called the excitation tensor. The components of G can be represented as a matrix, that 
in a local orthonormal frame (or Minkowski space-time with Cartesian coordinates) have values 
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(10) 



Then the identification 



G:=*F (11) 

is a linear map that takes ★F to G and provides a set of constitutive relations for the components of G in terms of 
those of F. Comparing Eq. ([TU]) with Eq. ([2]), it is clear that G = *F reduces to the trivial, vacuum, constitutive 
relations Ha — Ba and Da — Ea (where = — \ when c = 1). Thus we find that a trivial linear map recovers the 
correct constitutive relations in vacuum, but this will be extended to something non-trivial in the next section. 

We would like to stress that we take G to be a tensor, not a tensor density. Expressing the excitation as a 
tensor density is common in the literature, but using G = *F makes use of the metric structure and explicitly shows 
the space-time contributions encoded by which will be very useful when discussing materials and, in particular, 
transformation optics. 

Including an interaction term, the action is generalized to 

S = y^FAG + JAA. (12) 

In the interaction term, A is the connection 1-form, and J = Jaji-y is the charge-current 3-form. This action is 
invariant under a gauge transformation A — > A + d/ for some scalar function /. The 3-form J is related to the usual 
4- vector current j = j'^ by the volume dual, J ~ a;(j), or in component form 



(13) 
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What advantage is there to thinking in terms of the charge-current 3- form J rather than the 4-vector j? Again, the 
integral is only defined for forms of the same dimensionality as the manifold on which the integration is performed. 
The 3-form J may be integrated over a three dimensional hypersurface of the space-time manifold. Choosing a 
constant-time hypersurface is equivalent to a spatial 3-volume. Integrating J over this spatial 3-volume gives the 
charge enclosed. Integrating J over a 1-1-2 hypersurface corresponding to time and a spatial 2-surface gives the 
current flowing through the spatial surface. 

Equipped with the constraint equation dF = (the homogeneous Maxwell equations), and the constitutive Eq. 
(fTS)). varying the action of Eq. with respect to A gives 

dG = J. (14) 

This comprises the inhomogeneous Maxwell equations, although the reader may be more familiar with the expression 
obtained by taking the dual of both sides of Eq. (fT4|) . 



III. ELECTRODYNAMICS IN LINEAR MEDIA 



Electrodynamics in materials is somewhat more complicated than electrodynamics in vacuum. Here we expand 
on the brief introduction given in Ref. [l^ , whose main points are embedded here for a complete and self-contained 
discussion. The microscopic theory of electrodynamics would be a quantum field theory described by some complicated 
action that includes not only the electromagnetic fields, but also the various matter fields making up the material, 
along with their interactions and associated gauge fields [2^ [2§| . Considering the vast number of fields contained in 
a sample of ordinary material, it would be an impossible task to examine the full exact theory. Fortunately, in the 
thermodynamic limit, electrodynamics in media can be described by an effective theory. This generally comes in the 
form of a material dependent set of constitutive relations. 

The standard vector relations D = eE and H = n 
combined into an expression such as G^'^ ~ 
exploring electrodynamics in the absence of a metric [30|-|3^. However, since we assume the existence of a metric, we 
instead choose to retain the usual space-time notions of metric and Hodge dual *, and take the minimal approach of 
extending the trivial constitutive equation G = ^F in vacuum to a more general linear constitutive equation [l^ 



(where e and ^ ^ may be matrix- valued) are frequently 
2I, [mi. Such an expression can be useful, particularly when 



that in component form reads 



X, 



(*F)„^. 



(15) 



(16) 



The tensor x contains information on the dielectric material's properties, and can be thought of as representing an 
averaging over all the material contributions to an action that describes a more fundamental quantum field theory. 
The motivation for using this constitutive relation is to explicitly separate the space-time effects (i.e. the Plebanski 
relations) from the material effects, which will be useful for transformation optics. 

To retain the symmetry properties and usual notions of G and F, x must be independently antisymmetric on its 
first two and last two indices, and in vacuum x(*F) = Tk-F. Thus the classical vacuum is a perfect dielectric for which 
a trivial x describes the electrodynamics, and we extend this idea to a non-trivial x describing electrodynamics in 
arbitrary, linear, dielectric media. These conditions reduce the number of free parameters of x to 36. One may 
further decompose x into principle, skewon, and axion parts [26l [33l|. but we do not consider this here. Additional 
symmetry conditions may be imposed based on thermodynamic or energy conservation arguments, or by the lack 
of an observed directive effect in naturally occurring stationary materials [25l [3^ . Having recently entered an era 
of engineered materials which may incorporate active elements [35l [36| . however, we leave open the discussion of 
additional symmetries and consider the three conditions above to be the minimal requirements. The condition 
Xvac{*F) = ★F is sufficient to uniquely specify all components of x for the vacuum, they are 




(17) 
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Equation (|17l) expresses % as a matrix of matrices, the first two indices of X^^^ give the a/3 component of the 
large matrix, which is itself a matrix described by the second set of indices. The component values of x. for the 
vacuum are unique and independent of coordinate system. For a more general material, the component values can 
easily be determined by simply matching the results of the constitutive equation G = x(*F) with the usual flat-space 
constitutive relations in a particular coordinate system, as shown in Appendix]^ The components of the constitutive 
equation provide a set of six independent equations that can locally be collected in the form 

Ha = {fi-')a'B, + {ir),'E,, Da = {t)a'E, + {i2*)aB,, (18) 

where we use the notation d to denote a 3 x 3 matrix. Rearranging these to 

Ba^{il)a'H, + {i,)J'Et, Da^{e)J'E, + {i2)J'm. (19) 

gives the more familiar representation for the constitutive relations. These three-dimensional representations of the 
completely covariant Eq. (1151) are essentially equivalent, and it is a simple matter to switch between them using 

fi={fj-^^) \ e = - 72*M7i*, ii^~Mi*i 72=72*A- (20) 

One should be aware that these 3x3 matrices are not tensors, but simply components of x that have been collected 
into matrices. They could be made into tensors by incorporating the appropriate 3-dimensional Hodge dual of a 
space-like hypersurface. As they stand, Eqs. (fT5)l and (|19p are somewhat misleading, because while Ea and Ha are 
components of a 1-form, Da and are really selected components of the 2-forms Dabdx"" A dx^ and B^tdx'^ A dx''. 
Thus a more appropriate version of Eqs. pop would be something like 

{*^B)a = {il)a''Ht + {ii)a''Et, (21a) 

(*E^)a = {e)J'Et + {i2)a'Hb, (21b) 

where is the Hodge dual on the 3-dimensional space-like hypersurface. But this requires that we resolve the space- 
time into space and time components, selecting an observer to define a direction of time. The spatial hypersurface 
is then orthogonal to the selected direction of time. Transforming to the local frame of the selected observer we can 
make the identifications of Eq. (|A2[) and then give the constitutive equations a 3-dimensional representation. However, 
this requires greater care, and the inclusion of time transformations is not immediately evident. We will therefore 
continue with a manifestly 4-dimensional approach. 



IV. TRANSFORMATION OPTICS 

The completely covariant approach to transformation optics was outlined in [l6|, here we derive this method in 
greater detail. Start with an initial space-time M and field configuration {g,-k,x,F,G,3), where dF = 0, dG = J, 
and G = x(*F). The usual approach to transformation optics begins by imagining a coordinate transformation T 
that in some way "deforms" the manifold. For example, in the case of an electromagnetic cloak [1, [l^, we imagine a 
coordinate transformation that stretches out a hole in Minkowski space-time (with a point removed). 

However, certain subtleties involved in this picture should be clarified. The idea of stretching open a hole or 
otherwise deforming the space-time is non-physical; therefore it is more appropriate to imagine a map which takes 
M to an image M C M, as in Fig. [51 While it may be intuitive to think of the map T as given, the transformation 
of 2-forms (representing the fields of interest) requires instead a map T that describes how the image M is mapped 
to the original manifold. In certain situations it may be possible to set 7" = T~^, but in general T^^ may not even 
exist. Furthermore, while the transformation T acts to transform the fields, the space-time metric is transformed by 
a different map T, which here we take to be the identity. We will have more to say about this below. 

To make this more concrete, consider a pair of maps (T, T). We demand that % acts only on the metric via its 
puUback 

1: M ^ M (22) 
^*(g) = g. (23) 
On the other hand, we demand that the map T acts on the electromagnetic fields via its puUback 



r : M C M M 

r*(F) = F 



(24) 
(25) 
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(F,G,x) 
G = x{*F) 



(M,g,i) 

(F.G,x) 
'\ G = x(iF) 



FIG. 2: Under the map T the points of M are mapped to the image M. The electromagnetic fields are transformed by the 
puUback of T, T* ■ The metric, on the other hand, is transformed by the puUback of T, which for present purposes is the 
identity map. Thus M C M, and the pulled back fields, which exist on M, may be excluded from some region of M. 



as depicted in Fig. [2] Under (T, T), the metric and fields are transformed from the configuration (g, %, F, G, J), 
where dF = 0, dG = J, and G = x(*F) to a new configuration (g, x, F, G, J), where dF = 0, dG — J, and 
G-X(*F). 

The way to physically achieve such a transformation is to introduce some kind of material, just like introducing a 
dielectric between the plates of a parallel plate capacitor. Thus we know that the new configuration of electromagnetic 
fields must arise from a material distribution x f^nd must obey G = x(-a-F). The question of transformation optics is: 
given a transformation, is it possible to determine the x that will support the transformed fields? 

Consider Eqs. (|23l) and Eq. (|25l) in more detail. At a point a; G M we have 



(26) 



where the second line follows from the constitutive relation G = But from the constitutive relation G = 

we also have 



Gx^XxOixO T* (F 



nx)) 



(27) 



The right hand sides of Eqs. (|26p and (|27|) must be equal, but to clearly see what is going on we must consider the 
action of G^; on a bi- vector V^; e T^{M), so 



Xrix) ° *r(x) ° Fr(x} ° dTiVx) = [x T* (Frix))] (Vx) . 



(28) 



Letting A^^, be the Jacobian matrix of T (A^^ is the matrix representation of d7~) this may be written in component 
form as 



T{x 



r(x) 



T{x 



up 



T{x) 



(A%A'',) 



' X 1 



(29) 



where we have indicated explicitly where each tensor or Jacobian matrix is evaluated. Eliminating F and V^. from 
both sides and solving for x gives the final result 



r{x) 



T{x) 



T7] 



(30) 



In Eq. (|30p A^^ is the matrix inverse of A, both A and A^^ are evaluated at x, and in solving for x we have made 
use of the fact that on a 4-dimensional Lorentzian manifold, acting twice with ★ returns the negative, ★ ★ F = — F. 

What is ★ at the point x? The Hodge dual ★ is not directly transformed by the pullback of T. Rather, the metric 
is pulled back with T* and then * is computed from the pulled back metric. But notice that 



gx = 1*(gi(x)) = Si{x) o dT 



(31) 



depends on g at the point T(a;), which is not, in general, the same point as T{x). So x{x) depends not only on the 
point T{x), as explicitly shown in Eq. ([30]) . but also depends on the point T(x) through the evaluation of ★ at x. 
Equation ([50)1 is the core of transformation optics. Start with a given space-time with metric g and associated dual 
and with given dielectric material properties described by the tensor x- The initial space-time may be Minkowski 
and the initial dielectric may be the vacuum, but this is not necessary. We imagine a transformation that changes 
the fields in some way. We ask what x is required to physically achieve such a transformation. The answer is given 
by Eq. dSnH- 
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A few remarks are in order concerning Eq. (j30l) . For applications in transformation optics, T is generally taken to 
be the identity map, meaning that T(a;) = x and i = -k. However, this is not strictly necessary, and allowing % to be 
a more general map may be useful in the study of analog space-times [sH]. Furthermore, if in addition to T being the 
identity the initial space-time is vacuum, then since Xvac^ = *, 



(32) 



where the first ★ is evaluated at T(x) and everything else is evaluated at x [16]. Notice that the prescriptions of Eqs. 
([30|l or (f32|) are meaningful only for points x G M. So for transformations such as the electromagnetic cloak, where 
there is a hole in M, the material parameters inside the hole arc unspecified and completely arbitrary. In this way, 
any uncharged material may be hidden inside the cloak without affecting the behavior of the fields outside the cloak. 

One might be concerned by the step of assigning two different transformations, one to the metric and one to the fields, 
and argue that such a step is not allowed because it does not constitute a symmetry of the system. However, these 
are not maps between physically equivalent systems, clearly evidenced by the fact that one is vacuum and the other 
contains a material. Thus there is no symmetry that should be preserved. The point of transformation optics is that 
the transformed fields do not constitute a solution to Maxwell's equations in the original space-time but do constitute 
a solution to Maxwell's equations in the appropriate material. Thus the mappings involved are not mappings between 
equivalent solutions, but rather a method of generating new solutions from the original (untransformed) solution. 

Furthermore, one might be concerned that the puUback of the fields under T might not be contained in the image 
of the pullback of the metric under T. Ultimately, this map may only be defined locally. But since the final material 
parameters must ultimately be measured by a local observer with a local metric, the depiction in Fig. [5] can always 
be defined locally (neglecting irrelevant situations like transformation optics near singularities). 

Finally, we note that because of the symmetries of F and G, which have only 6 independent components each, 
and of X and 7*r, which have at most 36 independent components, we could re-express these results in terms of two 
6-dimensional vectors encoding the information of F and G and two 6x6 matrices encoding the information of x 
and *. This 6-dimensional representation is commonly encountered in the literature, and could be advantageous if 
performing calculations by hand because of the large number of terms that must be computed in Eq. pop . However, a 
modern computer algebra package can handle the full calculation with ease, and the manifestly 4-dimensional nature 
of the expressions is useful when further manipulations are performed. 



A. Examples 

The completely covariant approach to transformation optics developed above provides a concrete and powerful 
framework for analyzing any desired configuration of fields and linear dielectric materials in any space-time and with 
any relative velocities. Several previous examples have been given [T6| which demonstrate that the results obtained 
with the completely covariant approach agree with results obtained through other means. This section presents 
some additional examples to illustrate the applicability of the completely covariant approach to arbitrary dielectric 
velocities. 



Mixed Transformation 



Starting from vacuum Minkowski space-time, consider the transformation T defined by 

T{t' , x' , y' , z') — {t, x, y, z) — {ax't', x' , bx' + cy' , z') 



(33) 



The temporal part of this transformation has previously been considered in more detail p39| , but the spatial transfor- 
mation provides some additional complexity. As discussed above the fields are actually transformed by T rather than 
T, which we take to be 



T{t, X, y, z) = {t', x', y, z') ^ { —,x, - — z 



xj^O 



and the Jacobian matrix of T is 



A'',, 



ax 




ax^ 


'\ 





1 








_b 1 

c c 





^0 





V 



(34) 



(35) 
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Turning the crank on Eq. ([30]) . and comparing with Eq. (|A2[) . the components of x can be extracted and collected as 

62 + ^2-44 \^ j^-^^—i" -b 1 1 

n ° " 1 I acx \ ^2 / 

(36) 






Changing to the representation of Eq. (fT9|) by using Eq. pop results in 

4 /O 0\ 

(37) 

The time transformation is responsible for the appearance of the non-zero magneto-electric coupling terms [l6l, IseI] . 
It was shown in Ref. [16] that if e = /i is proportional to the identity matrix, then the magneto-electric coupling terms 
can be simply identified as a velocity, i.e. the magneto-electric coupling may arise solely from a material velocity. 
However, if e = /i is not proportional to the identity matrix, as in this example, then the magneto-electric coupling 
cannot always be interpreted as a simple material velocity. Thus in this example we must interpret these results as 
an anisotropic material with magneto-electric coupling, at rest with respect to the frame in which the fields are given. 
Indeed, it is not possible to boost to a uniformly translating frame in which the magneto-electric coupling vanishes, 
as we show below. Note that the time-dependent behavior of these material parameters make a validating calculation 
quite difficult. Such a calculation requires an analysis for both the initial and transformed mediums that is beyond 
the scope of this paper [s^. 

2. Moving Materials 

In a real world scenario, it is likely that the design parameters require the desired material to move with respect to 
the frame in which the field transformations are defined. In the Plebanski based approach outlined in Ref. |12| . the 
magneto-electric coupling term is interpreted as a velocity dictated by the transformation, rather than being a free 
parameter. As discussed in Ref. [T6j , a magneto-electric coupling can only be independently identified with a velocity 
if the material is isotropic. Therefore even if the material velocity can be arbitrarily tuned, it is not always possible 
to tune the velocity such that for a given transformation the material does not require inherent magneto-electric 
couplings. In the totally covariant method presented here, the material velocity is an arbitrary parameter compatible 
with any value of magneto-electric coupling. 

This section presents some examples of a transformation applied to a uniformly moving material in vacuum 
Minkowski space-time. For simplicity, suppose the resulting material moves in the a;-direction with some speed /3 
with respect to the laboratory frame in which the electromagnetic fields are known. A prime will denote an object 
in the material frame, while unprimed objects reside in the lab frame. A transformation to the material frame from 
the lab frame is done with a Lorentz boost such that for = , 



7 


-7/3 0\ 


7/3 


7 01 





10 





1/ 



L^u-\ 10' 
V 1/ 

where 7 = (l-/32)-i/2, xhe inverse boost, from the material frame to the lab frame, is obtained by setting /3 — /3, 
so that v'^ — L'^,,,v^ , where 



/7 7/3 ON 

7/3 7 

10 

V ly 



(39) 



An object with lowered indices boosts as o;^' ~ uj^L'^^, or = cj^'L^ . 
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The calculation now follows more or less straightforwardly from Eq. ([32]) . but includes a boost. Let x' be the material 
parameters in the rest frame of the material, the material parameters described in the lab frame are, schematically 

Xiab = (xO . (40) 



But the material parameters in the lab frame are precisely those obtained from Eq. (j32j), xiab — X- It follows that in 
the frame of the material, where the material parameters are to be measured 

Let the transformation again be that of Eq. p3p . Performing the calculation, extracting the components of e*, fi^^ , 
7i* and 72*, and converting to the representation of Eq. (fT9)) . results in 



a2x4(52^_c2)/32-c2(x + t/3)M 7^^ ':^^^^ 7^ ^ 

— 






o2a;4(62 + c2)/32 _c2(a; + t/3)" c^{x + tl3)(t + xP) - a'x'^ ^{h'' + c^) Oy 



(42) 



-Yl = -7,^ = - I 



where * indicates components that are obtained by the antisymmetry of the matrix. 

While these results are fairly complicated and perhaps do not offer much intuition regarding moving dielectrics, 
this example demonstrates the method by which material velocities may be easily incorporated into the theory, even 
for non-trivial transformations. Essentially, the calculation consists of calculating x in the frame in which the desired 
transformation is specified, and then transforming that result to the rest frame of the corresponding moving material. 
It can be readily seen that the limit /3 -> recovers the results of Eq. ([57)1 . 

As a somewhat simpler example, consider the transformation 

T{t,x,y,z) = {t' ,x' ,y' ,z') = (^at + bx,x, ^—J—'^^ ' (43) 



which is similar to that of Eq. p3p . but avoids the complication of time dependent parameters found in Eqs. p7p and 
(|42[) . For a stationary material, the equivalent material parameters are 

(44a) 




71 = 72^ = - 11. (44b) 



a 



,0 -1 0; 



It is tempting to associate the magneto-electric coupling with a simple material velocity in the x-direction, proportional 
to b/a. To check whether this is the case, suppose we take a material described by Eq. (|44ap for e = jj,, but with 
7i = 72 = 0, and moving in the x-direction with speed /3. In general, e = jl, and 71 = 72 will all now depend on /?. 
For low speeds, expanding to leading order in /3 finds e = jl unchanged to first order, but with 




71 = 72" - ^ I (a2d2 _ c2 - d2) I ^ (45) 

-iP) 

Thus interpreting the magneto-electric coupling of Eq. (I44b[) as a simple low-speed motion in the x-direction is 
incorrect, as this generates additional magneto-electric coupling components of the same order. However, 71^^ ~ 
cf3/a^cP indicates that including a y-component for the velocity might allow us to recover the desired results. Indeed 
this is the case, and it can be shown that a material with e = fi described by Eq. (|44al) and 71 = 72 = 0, moving with 
velocity components 



l-a2(i_c2)-/2(l_a2)' 
abc 

" l-a2(i_c2)-/2(l-a2)' 
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generates the results of Eqs. (|44|) to leading order in (3^ and fiy. Velocity induced corrections to e = (i are of second 
order, i.e. proportional to I3xl3y, and of course we are restricted to (^^ + ^2)i/2 ^ |^| ^< ^ 

3. Moving cloak 

Finally, consider a cloak moving in the x-direction with constant speed (3. What velocity induced corrections are 
required to achieve cloaking with respect to fields described in the laboratory frame? A square cloak is obtained from 
the transformation [s^l 

Tit, y, z) = it', A V\ = { , (47) 

for (si < X < S2), {—S2 < y < S2), and \y\ < \x\. From Eq. p2p . the well-known results for the equivalent material 
parameters of a stationary cloak are 

-^1* I (48) 



^r—^ 



X{S1-S2)'-' V 




with 71 = 72 = 0. Now assume that the cloak is moving. From Eq. (|4ip . it follows that in the frame of the cloak and 
to first order in the speed /3, the material parameters now have the additional magneto-electric couplings 

= ^ —J^1L_ I -a:^(2s2-s^)-\s^slj. (49) 

The moving cloak therefore requires inherent magneto-electric couplings designed to compensate for these velocity 
induced magneto-electric couplings, which are more complicated than might be expected from previous experience 
with isotropic materials. These corrections are likely to be small for ordinary non-relativistic velocities, but depend 
specifically on the other cloak parameters. As before, the corrections to e — fi are second order in /3. 

V. THE GORDON FORMALISM 

In a problem closely related to the methods of transformation optics, Gordon, in 1923, studied the possibility 
of representing electrodynamics in a material by an equivalent vacuum space-time i5;]. One motivation for such a 
representation is that it allows the material to be described by the action of Eq. p^ . with the constitutive relation 
G — ★F, where ★ is the Hodge dual for a vacuum space-time described by an effective, or "optical" metric. This is, in 
some sense, the "inverse" problem of transformation optics. In particular, suppose we are given a dielectric material 
described by x residing in a space-time with metric g and corresponding dual *, then from Eq. (1281) it may be possible 
to find ★ corresponding to an effective vacuum space-time, 

X* = Xvac* = *■ (50) 

To be more explicit, the initial setup consists of a material residing in a space-time, for which the relevant part of the 
action is 

S = I FAG^/ FAx(*F). (51) 

Similarly to Eq. ©, this is 



If this is to be regarded as a vacuum space-time with metric 7, then 



/ FAx(*F)=/ d^x^\g{F,-{^x*F)). (52) 
Jm Jm 



[ FAx(*F)== / FA*F= / d^x^\jiF,F). (53) 
Jm Jm Jm 
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The most general x has 36 independent components, but the metric only has 10 components, meaning that not 
every dielectric material can be represented by an equivalent space-time. To find the conditions that a material must 
satisfy to be representable as a vacuum space-time, return first to Eq. (|50|) . Calculating ★ for an arbitrary metric 
'y, and comparing with for arbitrary Xi it is readily observed that both e* and fj,~^ must be symmetric, while 
72* = —(71*)^ is traceless. This reduces the number of apparently independent equations to 21, plus one constraint. 
However, the equations are not linear and are given in a particular representation. Converting to the representation 
of Eq. p9p one finds that e = p., and that 72 = 'fi''' is not just traceless but antisymmetric. These, of course, are 
precisely the constraints expected from the Plebanski-De Felice equations. 

For a dielectric material to be representable by a curved vacuum space-time it must have e = fi, and 72 — 71^, 
which reduces the number of free parameters of x to 9. Unfortunately, this is insufficient to uniquely determine the 
10 free parameters of the metric. On a four dimensional Lorentzian manifold, every metric g is uniquely associated 
with a map * from 2-forms to 2-forms by 

*.r = ^^^a^.p5'^V^ (54) 

Notice that * is invariant to a scaling of the metric g ag, so given * it is only possible to determine the metric up 
to an overall scale factor. On the other hand this means that when trying to determine a particular 7 from a given *, 
one is free to choose a convenient representative from the equivalence class of conformally related space-times, such as 
the representative with 700 = —1. The situation may improve somewhat if there are sources present, J 0, since the 
Hodge dual of 1-forms and 3-forms is not invariant under conformal transformations of the metric. This may impose 
some additional constrains on the metric, but remains to be fully investigated. 

There exists another condition from which to obtain constraints on the material parameters. Recall that one of the 
defining characteristics of the Hodge dual on a 4-dimensional Lorentzian manifold is that **F = — F. It follows that 

(X*) (X*) = -*d2 (55) 

where we write id2 as the identity of the map from 2-forms to 2-forms. This does not provide any independent 
constraints on the material, but may serve as a useful check when performing calculations. 

It is not our purpose to find solutions for a particular distribution of materials, and it is sufficient to remark that, for 
the case of an isotropic material with vanishing magneto-electric coupling, it is relatively easy to check that Gordon's 
optical metric 

l^..=9^.. + {l-e^l)Uf,u,, = g""" - {e - ^) u" , (56) 

where is the 4-velocity of the material, agrees with our method. For an observer at rest with the material, 
= (1, 0, 0, 0), while for a material in relative motion to the observer, the x in Eq. (|50l) is first transformed with the 
appropriate boost. 

VI. CONCLUSIONS 

We have continued the development of a fully covariant and manifestly 4-dimensional formalism for transformation 
optics in general linear materials that began in Ref. (1(|. The main result of this construction, Eq. (j30[) . represents 
a sort of master equation for transformation optics of general, linear, nondispersive media. The benefit of this 
approach is that it is valid for arbitrary background space-times jioj , arbitrary initial material [Tsj , and for arbitrary 
transformations. The motivation for this construction was the observation that the Plebanski-De Felice equations that 
are frequently used in transformation optics and identified with material velocities are only valid in the low velocity 
limit for isotropic materials. 

To arrive at the completely covariant method it was first necessary to find a modest generalization of electrodynamics 
in linear dielectric media. This came by way of a modification to the constitutive relations to the form G = x(*F). 
Within these relations the vacuum is identified as being a perfect linear dielectric with relative permittivity and 
permeability equal to 1 - corresponding to a trivial x ^ ^-nd is therefore treated just like any other dielectric media. 

Next, the interpretation of transformation optics was clarified by realizing that the system is being transformed 
by two distinct maps: one applied to the metric and a different one applied to the electromagnetic fields; because a 
transformation of the fields is physically realized by changing the dielectric material properties (e.g. by introducing a 
non-vacuum dielectric material to a vacuum region) but generally not the metric. 

As an illustration of the completely covariant approach we have examined the results of a transformation when the 
resulting material is constrained to move with arbitrary uniform velocity - a physically realistic scenario that has 
not, as far as we know, been previously considered. This is a departure from previous results, where the material 
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velocity has been dictated by the transformation or has been constrained to the nonrelativistic hmit. It is furthermore 
shown that the magneto-electric coupling does not arise as straightforwardly from a material velocity as expected from 
previous interpretations of such couplings. For the particular transformation considered it was found that it is in fact 
possible to approximately get magneto-electric couplings from a material velocity, but only in the low velocity limit 
where corrections to the permeability and permittivity are neglected, and even then the velocity was in an unexpected 
direction. It was shown that for a moving cloak, magneto-electric couplings on the order of the velocity have to be 
designed into the cloak to compensate for the velocity-induced couplings. 

It was shown in Sec. IVlthat our formalism generalizes, in a covariant way, Gordon's optical metric Q, which turns 
out to be essentially the inverse problem of transformation optics. Namely, given a dielectric, can we find an equivalent 
vacuum space-time? The answer appears to be yes, but only if the dielectric satisfies certain properties, and even 
then only up to a conformal factor of the metric. 

It is worth reiterating that the component form of the constitutive relations in the covariant formalism are of the 
form of Eq. (fT8|) rather than Eq. (|T9l) . This leads to a slightly different meaning for the 3x3 matrices e*, p,~^, 71*, 
and 72* than what would be expected from the relations Eq. (|19p . One may readily switch back and forth between 
the two representations for the constitutive relations via Eq. (|20p . Keep in mind, however, that the 3x3 matrices of 



Eq. (jl8p are merely selected components of the main quantity of interest, the covariant („) -tensor x- 
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Appendix A: Cartesian Coordinates 



In Cartesian coordinates, the field strength tensor is defined as 



F = Exdx Adt + Eydy A dt + E^dz Adt + B^dy Adz- Bydx Adz + B^dx A dy. 



(Al) 



This may be represented in matrix form by Eq. ([2]). To recover the usual component relations of Eq. (1181) the 
constitutive equation G = x(*F) allows us to make the identification [iGj 



ap _ 














-72x 















(A2) 



where the * indicates entries that arc antisymmetric on either the first or second set of indices on x^^^- 

In cylindrical coordinates (i, r, 0^ z), the Minkowski metric tensor is = diag{—l, 1, r^, 1). Because we have c — 1, 
time and space are measured in the same units; meaning that E" and i?" are both measured in units of V/m and 
the units for F are V ■ m. Changing to cylindrical coordinates does not change the units of F, so 



E^rdr Adt + rEodO Adt + EAz Adt + rBrdO Adz- Badr A dz + rBAr A dO. 



(A3) 



This definition for the components of F in cylindrical coordinates assumes that all E"' and i?" are still measured in 
units of V/m, despite the dimensionless coordinate 9. Similar arguments hold for G. These are written in matrix 
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form as 




and Gaiy = 








rHg, 
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-rHg 


-rD, 







Dg 


~rDr 




(A4) 



Matching the constitutive equation G = x(*F) for Minkowski space-time in cyhndrical coordinates resuhs in the 
identification 
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(A5) 



The form of the components in both of Eqs. (jA4l) and (lASp could equally well be determined by transforming the 
Cartesian versions, Eqs. (jlOp . and (jA2p . with the appropriate Cartesian-cylindrical transformation matrix. 

In spherical coordinates {t,r,9,ip), the Minkowski metric tensor has components 17^1, — diag{— 1,1, r^,r^ siri^ 6). 
This time there are two dimensionless coordinates, so 



F = Erdr Adt + rEgdO A dt + r sin 9E^dip A dt + sin OBrdO A dip ~ r sin OBgdr A dip + rB^dr A d9 
(every d^ gets a factor of r, every dip gets a factor of r sin^?), which is written in matrix form as 

/ -Er -rEg, 



F„ 
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(A6) 



(A7) 



Letting s 
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sm9, the same identification procedure as before leads to 
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(A8) 



J 



Like the cylindrical result, the spherical result could be obtained equally well by transforming the Cartesian Eqs. 
pop , and (|A2p , with the appropriate Cartesian-spherical transformation matrix. 
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Appendix B: Fullback and Pushforward 



Suppose two manifolds are related by a smooth map (f : M ^ N. Thus (p maps the points of M to an image in N, 
p E M ^ (p{p) G N. Given a function / on N, the pullhack of / by 95, denoted (p* f is the composite function f o tp. 
So the pullback (f* takes a function / on iV to a (p-related function (f*f on M. 

The differential of (/? at a point p & M defines a map of tangent spaces dtp : TpM — >■ T^(^p)N. Given a vector 
Vp G TpM, the pushforward of Vp assigns a (^-related vector dipiVp) G T^^p-^N, such that 



MvUp)f = Vp{<p*f) 

for some / on N. In particular, given coordinate charts {x^} on M and {y^} on A'' this is 

(My,K^p,^f = Vp^-^{^*f) 



dip° 



A;. 

9y« 



Thus the pushforward of Vp is 



where 



i.MVv))l(p) = f^\Vp^, 



(Bl) 

(B2) 
(B3) 
(B4) 

(B5) 
(B6) 



is the Jacobian matrix of the map p evaluated at p. Note that the pushforward of a vector field is not always defined. 
For example, if p is not injective then for two distinct points p,q € M such that p{p) = p{q) e N, the vectors Vp and 
Vq would be pushed to the same point in N, but there is no reason why d(p{Vp) should equal dp>{Vq). However, if tp 
is a diffeomorphism, then the pushforward of a vector field is well-defined. 

Like a function (which is a 0-form), a 1-form on N can be pulled back to a 1-form on M such that 



{p*w) {Vp)=u;{dp{Vp)). 
In particular, expressing this in coordinate charts as above, we find 



da;' 



Which simplifies to 



((^*W^(p)) I 



dpP 




or 



{p*oj)a =a;/3A^„, 
where A^^ is again the Jacobian matrix of ip evaluated at p. 
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(B8) 



(B9) 
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